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ABSTRACT: Among the inflationary models based on string theory, the D3/D7 model has the 
advantage that the flatness of the inflaton potential can be protected even with moduli stabilization. 
However, the Abrikosov-Nielsen-Olesen BPS cosmic strings produced at the end of original D3/D7 
inflation lead to an additional contribution to the CMB anisotropy. To make this contribution 
consistent with the WMAP results one needs an extremely small gauge coupling in the effective 
D-term inflation model. Such couplings may be difficult to justify in string theory. Here we 
develop a generalized version of the D3/D7 brane model, which leads to semilocal strings, instead 
of the topologically stable ANO cosmic strings. We show that the semilocal strings have unbroken 
supersymmetry when embedded into supergravity with Fl terms. The energy of these strings 
is independent of their thickness. We confirm the existing arguments that strings of such type 
disappear soon after their formation and do not pose any cosmological problems, for any value of 
the gauge coupling. This should simplify the task of constructing fully realistic models of D3/D7 
inflation. 
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1. Introduction 

Despite the recent period of intense research activity on inflationary models in string theory, 
it is still a challenging task to derive a realistic model of inflation consistent with the 
stabilization of all string moduli. A particularly important problem is to reconcile the 
original idea of brane inflation jl], ||] with an appropriate mechanism of volume stabilization 
so as to make it a valid model of string theory. One of the most developed models in this 
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context was suggested in Q based on the mechanisms of volume stabilization proposed 
in H). Recently a particular class of models in string theory with moduli stabilization 
was constructed in Q. A major problem associated with inflationary models in warped 
compactifications of type IIB string theory has to do with the fact that the inflaton in 
these models (the position of the D3 brane deep inside the warped throat geometry) is 
a conformally coupled scalar in the effective four-dimensional near de Sitter geometry. 
Therefore, the mass of such a scalar, rrvi, is close to 2H 2 = ^R, where H 2 is the Hubble 
parameter and R is the curvature of de Sitter space. This does not meet the observational 
requirement wA ~ 10~ 2 H 2 , as was uncovered in |J, where some possible resolutions of 
this problem were suggested. Additional attempts to solve this problem in closely related 
settings were recently proposed in Q. New ideas for string cosmology based on non-linear 
Dirac-Born-Infeld brane actions were developed in ||. 

On the other hand, it has been found recently || [lOL 11, 12] that a particular brane 



inflation model, the D3/D7 model Q, features a shift symmetry for the inflaton field. The 
inflaton in this model is related to the position of a D3 brane with respect to a D7 brane 
in a K3 x ^- compactification of type IIB string theory, where Z2 involves both orbifold as 
well as orientifold operations, Z2 = I2 • ^ • (— 1)^ L . 1 The inflaton shift symmetry follows 



from the special geometry in such models, as described in ref. [13]. 



The D3/D7 model H has an effective description as a hybrid inflation [14], or, more 



specifically, as a D-term inflation model [15, 16]. The flat direction of the inflaton potential 



is associated with the shift symmetry, which is preserved after the stabilization of the 
volume modulus. It is broken spontaneously only at the quantum level. This symmetry 
protects the small inflaton mass during inflation. This is an important advantage, since in 
other models of inflation in string theory it is difficult to protect the inflaton field in the 
early universe from acquiring an unacceptably large mass. 

In its original form, however, the D3/D7 model with U(l) gauge symmetry has 
another problem. At the end of inflation, cosmic strings form, giving a contribution to 



the primordial density perturbations [17|. Such stringy perturbations do not lead to the 
multiple peak structure in the spectrum of the CMB anisotropies. The data from WMAP 
show that the contribution of the cosmic strings to the CMB anisotropy should not exceed 
a few percent of the standard inflationary perturbations [jl^, [l{J. In a general class of 
brane inflation models, the production, spectrum and the evolution of cosmic strings were 



studied in f20j, |2l], [2g, |j| ||, |T§. 



It is interesting to note that the model (3| of D3/D3 inflation in the highly warped 
throat has a natural mechanism for suppressing the contribution of the cosmic strings to 
the CMB anisotropy ||: The effective tension of the produced strings and, therefore, their 
contribution to metric perturbations, is exponentially suppressed due to the large warp 
factor. The production of very light cosmic strings in this model and in its generalizations 
was recently studied in |25|] . Note that in the warped region of the geometry there is an 



1 We use the standard notation I2 and f2 to denote, respectively, the orbifold and orientifold action, 
whereas (— 1) Fl reverses the left-moving spacetime fermion number. 
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m 2 ~ H 2 problem with inflation. On the other hand, large warping is known to relieve 
the problem with cosmic strings. In the D3/D7 model, by contrast, we do not have the 
m 2 ~ H 2 problem, but we also do not have an automatic solution of the cosmic string 
problem. 

There are many ways to address this problem. One may change the theory in such a 
way as to make the potential somewhat asymmetric, so that stable strings are not formed 
1 26 , 27]. Alternatively, one can consider models with extremely small coupling constants, 
which suppresses the stringy contribution [p8| , 17, ^9|. In particular it was shown in [17] 
that if one assumes, following IS], that the contribution of cosmic strings to the amplitude 
of metric perturbations should be smaller than O(10~ 2 ) of the inflationary contribution, 
then the effective g 2 must be smaller than 10 -10 . Meanwhile if one follows recent work [19] 
and assumes that the stringy contribution should not exceed O(10 -1 ) of the inflationary 
contribution, then the constraint on g 2 which can be obtained following Ref. [17| becomes 
somewhat less stringent, but it still remains quite strong: g 2 <^ 10~ 7 . In the context of the 
D3/D7 model with volume stabilization, one expects that the effective g 2 comes from the 
D3 brane gauge coupling. It may then be rather difficult to get g 2 <^ 10~ 7 in this context. 

A simple and elegant method for solving the cosmic string problem in D-term inflation 
models has recently been proposed in [3C] and 16]. If one simply replicates the charged 
matter multiplets responsible for the string formation so as to obtain a non-trivial global 
symmetry such as SU(2), one can render the vacuum manifold simply connected, and the 



relevant cosmic strings become so-called semilocal strings [31]-|^6[. The stability proper- 
ties of semilocal strings cannot be inferred from topological considerations alone, but also 
depend on the values of the couplings in the underlying Lagrangian |31|- |]35|] . Numerical 
simulations of their evolution suggest, that they do not pose any cosmological problems 
1 37]. As an alternative way to circumvent the cosmic string problem, it was suggested in 
1 16] to embed the D-term inflation model into some model with larger gauge symmetry so 
that the strings would not be topologically stable. 

Following these lines, we will, in this paper, generalize the original version of the 
D3/D7 inflationary model Q in such a way as to include a larger set of global and/or local 
symmetries. In order to do so, we first go back to the original model Q and discuss the 
precise relation between the brane construction and the associated four dimensional N = 2 
supersymmetric gauge theory. In this discussion, we will also carefully take into account 
the volume modulus dependence at each step so as to update the model with regard to 
the more recent work on stabilization. The purpose of revisiting the original model |2| is 
to show how exactly the surviving interacting U{\) gauge group of the effective D-term 
inflation model is related to the two £7(l)'s from the D3 and D7 worldvolume theories. 
This will give us a much better understanding of the original model and will allow us to 
proceed with generalized models with higher symmetries resulting from brane systems with 
multiple branes. Furthermore, when we embed this discussion into the compact model of 
Q but now with stacks of coincident D3 or D7 branes, we will show that the underlying 
F-theory picture correctly predicts the local and global symmetries one expects for the 
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cosmology of the model. 

Fluxes and non-perturbative effects, providing stabilization of moduli in string theory 
1 38 , U , as well as embedding of the relevant 4-dimensional gauge theory with FI term into 
4-dimensional supergravity break M = 2 supersymmetry to J\f = 1 supersymmetry. 
The embedding of the standard Abrikosov-Nielsen-Olesen (ANO) cosmic strings into M = 1 
supergravity with FI terms have been studied recently in |23|. They were shown to saturate 



the BPS bound and to preserve one half of the original N = 1 local supersymmetry [23]. 
Their possible relation to Dl strings and wrapped D3 branes in string theory was discussed 
in H H,|T|. 

As we show in this paper, the semilocal strings relevant for our generalized D3/D7 
inflationary model with global symmetries can likewise be realized as 1 /2-BPS solutions in 
M = 1 supergravity. More precisely, we show that there are families of 1/2-supersymmetric 
semilocal string solutions parameterized by an undetermined integration constant that 
controls the thickness (or the "spreading" ) of the flux tube in the string core. These findings 
are in agreement with earlier results in related non-supersymmetric models [32, 3^, 34] and 



imply that it does not cost any energy for the semilocal strings of our model to change 
their thickness. 

Finally, we will discuss the production and the cosmological evolution of semilocal 
strings after inflation in our scenario. Unlike the topologically stable ANO strings, which 
form loops or can be infinitely long, the semilocal strings formed after inflation look like a 
collection of loops and open string segments of finite length with monopoles at their ends. 
This fact, together with the possibility of their spreading at no energy cost, may lead to the 
absence of long semilocal strings, and, as a result, to the absence of the string-related large- 
scale metric perturbations. Indeed, numerical simulations suggest that semilocal strings 
disappear soon after their formation |j37f| . However, these simulations were performed 
in the context of a slightly different model, without taking into account the cosmological 
evolution. For a full quantum mechanical investigation of this issue one can use the methods 



of lattice simulations of symmetry breaking after hybrid inflation [39|. We hope to perform 



these calculations in a separate publication. In this paper we restrict ourselves to a simple 



qualitative analysis of the situation. We will argue that the main conclusions of Ref. [37] 



should be valid for the semilocal strings produced after inflation in the model studied in 



|3C] and in our paper. This suggests that the generalized D3/D7 model developed in our 
paper does not have any problems with cosmological perturbations produced by strings, 
and therefore there is no need to require that the coupling constant g 2 must be very small. 

Thus, the purpose of this paper is threefold: i) To develop the generalized D3/D7 
brane model including the volume modulus as well as some global symmetries in addition 
to the local U(l) gauge symmetry; such models will have semilocal cosmic strings instead 
of ANO strings, ii) To establish the unbroken supersymmetry of these semilocal strings 
embedded into M = 1 supergravity with FI terms for arbitrary "spread" parameter, iii) To 
critically analyse and confirm the arguments existing in the literature |k]] that the models 
with semilocal strings have no cosmological problems. 
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The organization of this paper is as follows: In section 2 we take a single D3 and a 
single D7 that is wrapped on a K3 manifold with non-compact orthogonal directions. We 
study the precise world volume actions of the branes when the system is kept in a curved 
background with appropriate fluxes turned on. This analysis will serve as a dictionary 
between the world volume coordinates and the known P-term action studied earlier in |17J . 
In section 2 we do not take the back-reactions of branes and fluxes on the geometry into 
account. This is addressed in section 3 where we also consider a fully compact scenario 
using the underlying F-theory picture that was developed in Q. We show that the F- 
theory curve correctly reproduces the local and the global symmetries for the system when 
we have multiple D3 and D7 branes. The back-reactions of branes and fluxes are shown to 
be given by warp factors that are typically of order one when the K3 manifold is fixed at a 
large volume. In section 4 we come to the analysis of semi-local strings in this model with 
U(l) gauge symmetry and a global SU(2) symmetry. Earlier studies of semi-local strings 
in the literature did not address the precise connection with local M = 1 supersymmetry. 
We show that the semi- local strings can arise as 1 /2 BPS solitons in the extended Abelian 
Higgs model embedded in J\f = 1 supergravity. Finally in section 5 we discuss the cosmology 
of the model. We discuss the issue of formation of semi-local strings and their effect on 
the large scale density perturbations of the universe. We argue that these strings do not 
cause any problems with cosmological perturbations at any coupling. Therefore, we are 
not restricted to keep the coupling constant g 2 very small. This may help to construct 
fully realistic models of D3/D7 inflation. 



2. The D3/D7 model, P-term model and the volume modulus 

In the first version of the D3/D7 model of inflation [Q] the internal manifold was assumed to 
be compact, but the mechanism of stabilization of the volume modulus was not discussed. 
However, later in |^, |12[ the volume modulus and its stabilization in the "inflaton 
trench" were included into the effective D-term inflation model. Here we will study the 
dependence on the volume modulus in the D3/D7 brane system by carefully keeping track 
of the powers of the volume, and we will also clarify the origin of the effective 4-dimensional 
model. Postponing the discussion of six compact dimensions to section 3, we will first 
consider putting the D3/D7 system on a compact K3 with the D7 wrapped on the K3. 
We will also give an explicit map between the usual D-brane variables and the variables 
used in cosmology from d = 4 supergravity as in [|l7]] . 

The basic features of the D3/D7 inflationary model ||] can be easily understood from 
the general properties of Dp/ D(p + 4) brane systems (see e.g. ^(|). As worked out in [40] 



for the special case p = 5, 5 — 9 strings have charge +1 with respect to the U(l) gauge 
field on the 9-brane and charge —1 with respect to the U(l) gauge field on the 5-brane, 
whereas for 9 — 5 strings the signs are reversed due to their opposite orientation. Two 
T-dualities then yield the system with p = 3 which we are interested in. We thus have a 
D3/D7 system with 3 — 7 and 7 — 3 strings stretched between the D3 and D7 branes with 
completely analogous charge assignments. 
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In the case when the fluxes of the D7 gauge field in the directions orthogonal to the D3 
brane are self-dual, the masses of the 3 — 7 and 7 — 3 strings depend only on the distance 
between the D3 and D7 branes, and the system has N = 2 unbroken supersymmetry. In the 
case when the D7 fluxes are not self-dual, supersymmetry is spontaneously broken and the 
masses of the stretched strings have an additional contribution proportional to the amount 
by which the world volume fluxes are not self dual {{6\ — 62 ) 2 in the notation of P])- The 
sign of this contribution depends on the charges of the strings: one string gets a negative 
contribution and the other one gets a positive contribution, in addition to the positive 
mass coming from the distance between branes. As a result, the positively charged string 
may become tachyonic at the critical point where the distance between the branes precisely 
cancels the contribution due to the non-self-dual flux on the D7 brane. This allows one 
to identify the D3/D7 brane system with the M = 2 version [28, 17 1 of the gauge theory 



of the D-term inflation model [15], with one hypermultiplet [|[ and a Fayet-Iliopoulos (FI) 
term. 

The splitting of mass in the hypermultiplet leads to the logarithmic correction to 
the gauge theory potential, which can also be reproduced by the calculation of the one- 
loop amplitude for the effective interaction between the branes. This field theory one-loop 
correction will have an additional factor in case that one has more than one hypermultiplet. 



The end of the waterfall stage of standard D-term hybrid inflation [14] solves the D- 
flatness condition, and one ends up with an unbroken supersymmetry and a spontaneously 
broken U(l) gauge symmetry. In the D3/D7 brane system, this final state corresponds to 
a bound state, in which the D3 brane is dissolved as a deformed Abelian instanton on the 
U(l) DBI theory on the D7-brane such that the Chern-Simons coupling on the D7 brane 
is 2 

Scs = ~^ [ C±I\TI\T=\ C 4 . (2.1) 
16?H J D7 J D3 

Here, T mn = (cL4 — B) mn has a contribution from both the D7 worldvolume gauge field 
and the background NSNS two- form (in the directions orthogonal to the D3). The NSNS 
two- form plays the role of a noncommutative deformation parameter in this context. The 
anti-self-dual part of the deformation parameter is related to the constant part of (J 7 — *J r ) 
on the D7-brane. Only when the deformation parameter is non-vanishing does the Abelian 



instanton solution have finite energy [41|. The D3/D7 bound state has an unbroken super- 
symmetry, derived from At-symmetry on the D7-brane. The endpoint vacuum is therefore 
described by a non-marginal bound state of D3 and D7-branes that corresponds to the 
Higgs phase of the gauge theory of the D3/D7 system with the FI term provided by the 
nonvanishing 9\ —62- 

The low energy effective action (i.e. to leading order in g s and a') for a single D3 brane 



For the compact case, i.e. when we have K3xT 2 /Z 2 , we also need to incorporate the contribution of the 
A-roof genus terms in Scs- In that case, the instanton number f KS T A T should be ( 16 M j K3 tr TZ A 
TZ + 1), where TZ is the curvature two form. 
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and a single D7 brane in flat space is 



S 



d A x 



+ / d 4 x 



1 



ST 



2gZ(2iT(x 



- 2 {(d,X*) 2 + (d m X*f) 



2^(2™') 
2ira' 



- {{d^f + (d,Y' m ) 2 ) 



\X\ 



93 



(x^ A x) 2 



+ Scs ■ 



(2.2) 



The indices are as follows: [i runs from through 3; i denotes directions 4 and 5; m extends 
from 6 to 9; and S, T will include both the \i and m directions. The primed (unprimed) 
fields are the light open string degrees of freedom on the 3-brane (7-brane). The x doublet 
arises as the lightest degree of freedom from stretched strings between the D3-brane and the 
parallel D7-brane. That x is generically heavy is indicated by its mass which is proportional 
to {X % — X' % ) 2 . The covariant derivative is D^x = + iAp — iA'^)x- Note that x has 
units L^ 1 appropriate for a scalar in four dimensions. Note also that the fields X 1 have 
units of L, and has units of L _1 in all dimensions. Recall that g 2 = (2ir) p ~ 2 g s a'( p ~ 3 ^ 2 . 
The last term contains the various Chern-Simons terms on the two branes. 

We would like to place this system in a background that has some B mn field, a Ramond- 
Ramond F§ = dC^ turned on (consistent with the underlying orbifold as well as orientifold 
operation) and has the curved metric 3 



ds 2 = R- Q {x)gf v dx p dx v + R 2 {x){ds\ 



7: 



? 2 + ds 2 K %) 



(2.3) 



Here, g^(x) is the four-dimensional metric in the Einstein frame. The function R(x) is 
a function of the noncompact directions x p and is to represent the overall length scale of 
the compact K3 4 . The powers of R(x) in each part of the metric are chosen so that the 
four dimensional gravitational action is in Einstein frame ||. The function R(x) should be 
thought of as a very weak function of x (we will drop its derivatives) which we would aim 
to stabilize by wrapping branes on some cycles in the fully compact K3 x T 2 /Z2, invoking 
the mechanisms of 0]. 

2.1 The D7 action 



First, we will consider the D7 part of the action. In the curved background ([^), the first 
line of ( [2.2j ) modifies to 



S D7 = J d 4 xd 4 y^g^R- 12 ^gl^R 4 



1 



- -7-2 (FfivFpag 1 ^ ' g^E + F mn^F rs#X3#jO-^ 4 ) 

^5V 



2g 2 .(2ira') 2 E 13 



+ H7 



[2-nol 
27 



C( 4 ) A T A T , 



(2.4) 



3 In this section, we are only compactifying on the K3 directions so as to avoid having to cancel the 
charges of the 7-branes and 3-branes. Note that we are also ignoring the warp factors arising from the 



back-reactions of branes and fluxes on the geometry. In the next section we will consider a fully compact 
ario. 

Here, for simplicity, we keep track of only one of the several Kahler moduli of the compact space. 



scenario 

4 
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where we have neglected the fluctuations of Xi in the K3 directions and exhibited the 
powers of R(x) explicitly 5 . The covariant 2-form on the D7 brane consists of 2 parts 



T = F-B = dA-B. (2.5) 

Here, F mn is the field strength of the vector field A m which lives on the brane and B mn is 
the pullback of the space-time NSNS two-form field to the worldvolume of the D7 brane. 
We have also added the relevant Chern-Simons piece induced by the background RR field. 
Now we take J d^y^JgKz = Vk3 to be the volume of some fixed K3. Carrying out the 
integral over the K3 then yields 



f-l\ /' ] 



i\2 



A' 3 



where the coupling constants are 



—g- = T 7 (2,a ) VksR = = g • ^ 

Note that is part of the background, not from the D3. The T factors in this expression 
are along K3, and Cu\ is extended in the flat, non-compact directions. The background 
four-form is Cu\ = g s 7rR~ 12 Vol^/2. The Chern-Simons term then combines with the 
kinetic term for the gauge field in the K3 directions to make the second line 

-12 /"I 



Vol R-^l \ / jf-a*^- , (2.8) 

\°97 / JK3 

where (J- — *J-) = J-~ . Here, the * is the Hodge star on the fixed volume K3. One 
may think that would not point in the non-compact directions, but recalling that the 
SL{2, r L) invariant five form needs to be self (Hodge) dual in ten dimensions tells us that 
there's C( 4 ) in all ten directions. Either way, the action is now 



S D7 = I d'x^ ( ~^(F, u f - ,J _„ 2 (d» Xn2 

1 



Vol {A) R~ 12 ^ f J 7 ' A *T~ . (2.9) 



J7 JK3 



If one wanted to add more coincident D7s at this point, this expression would be suitably 
generalized by making the connection An a U (Nj) gauge field, letting the field X 1 be U{Nj) 
adjoint valued, giving it the appropriate covariant derivative and adding an overall trace 
to the entire expression. 



5 As mentioned earlier, there is an additional coupling coming from the tr 1Z A 1Z term on the seven 
branes. This coupling comes with a relative minus sign. We will henceforth assume that we have chosen the 
instanton numbers in such a way as to cancel this contribution. As a bonus, this will also avoid the issues 
of enhangons etc. jl2| in the system. For the compact case, which we discuss later, these complications will 
not be there. 
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2.2 The D3 and stretched strings 

Now let's put the D3 on the same metric and RR background. The expression becomes 



Sd3 = / d A xyJ-gE 



4g 2 



(2.10) 



~±( F > )2 R 

As far as the hypermultiplet is concerned, there are now three modifications with respect 



(F - )2 -S(2W)5(( 8 - Y ' ,)2 + (i ''- i ''"' )2 ) 



to the simple flat space action ( |2.2| ): 

(i) The mass for the hypers (from (X 1 — X' 1 ) 2 ) picks up some factors of R. 

(ii) The g| that appears in front of the D-term changes since the hypers are also charged 
under the gauge field on the D7 6 . 

(iii) There is also an overall factor of R~ 12 . 
Putting everything together, the new action is 



1 

471 



" A „2 C^A">) 



4^3 

R- 4 



2^(2™' 



/7 ^-^3 



R 



((d.x'r + (d,Y ,m f) 



2^(2vra' 



-10 



2vra' 



IXI 



■{x^ A x? 



(2-11) 



In this equation, the hypermultiplet covariant derivative is still D^x = (dp + iA^ — iA')x- 
We perform a gauge field redefinition 



gW, = A fl -A' fX gW = pA„ + — A' g 2 = g 2 + ~g 2 

93 93 

under which -^F 2 + \F' 2 = F'h, + The scalars must also be rotated, 

^((|^ + |x- ) + i (|^ + |x- 



(2-12) 



gS' 



(2-13) 



after which their kinetic term has the standard form. In order to also bring x m ^o the 
standard form, we must rescale, x ~~ R^X- The action at this point is 



S = J d A xy/-g E 



-\{F w ?-\{F w ,) 2 -\dS\ 2 -\dS>\ 2 



R 



-12 



1 

8gj Jk3 



(2-14) 



where we have ignored the motion of the D3 brane parallel to the wrapped direction 
of the D7 brane. We have also neglected derivatives of R(x). With these gauge fields, 
DaX = (Pu + W^u)Xi an d all fields have the standard normalization of M = 1 chiral 
super fields. 



'Note that this contribution is not present when the K3 is of infinite volume, as in (2.2). 
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2.3 P-term action 



Now since this system is known to have rigid N = 2 supersymmetry (broken to M = 1 
when coupled to gravity) in d = 4, there must be some way to express this action in terms 
of standard J\f = 1 language. The suggested expression in terms of (rigid) J\f = 1 chiral 



superfields is the so called P-term model [17| 

C = ~h% - \D^ + \ 2 - \D^^ 2 -2g 2 ( |.S'|-(|<T» + |- ! |<I>_|-) + 



$ $ - it 

+ 2 



-^d<i> + i 2 -i<i>_i 2 -e3) 2 . (2.i5) 

In this expression, the fields are as follows. The complex fields and sit in two chiral 
multiplets of opposite charge with covariant derivatives D^± = (d tl ±igW lJ )<&±. The field 
S is neutral under the gauged U(l) and likewise sits in an N = 1 chiral multiplet. In N = 2 
language, the pair ($+, $_) should be viewed as an M=2 hypermultiplet charged under a 
U{\) gauged by the M = 2 vector multiplet (W^, 5). <3?+ and form a doublet under the 
M = 2 R-symmetry group SU (2)r. In J\f = 1 language, the potential can be obtained from 
the superpotential = \/255($ + $_ - £ + /2) and D-term D = |$ + | 2 - |<I>_| 2 - £ 3 . The 
coupling constant g in the superpotential is the same as the U(l) gauge coupling so that 
the system is in its N=2 limit. The U(l) Fayet-Iliopoulos term ^3 has been included. Note 
that the fields in the P-term expression are chiral superfields with the standard kinetic 
terms of the type D^&D 11 ®*. 

2.4 Comparing the actions 



The action in eq. ( 2.14| ) is to be compared with the P-term action, eq. ( 2.15 ). Since the 



fields in the brane expression are already normalized as N = 1 chiral fields, the comparison 
is straightforward. The neutral scalar S is to be identified in the two expressions. The 
hypers should be identified as xi = ^+ an d X2 = so that the covariant derivative 
behaves correctly, 

l^xl 2 = + igW^xl 2 -^\(d^ + igW^ + \ 2 + Kfy - igW^.l 2 . (2.16) 

Just as in the P-term action in section 2.3, the fields (W^,S) should be thought of as an 
N=2 U(l) vector multiplet. The other gauge multiplet, (W' S') has completely decoupled 
and is irrelevant in the P-term model. From the brane picture, it is the decoupled center 
of mass degree of freedom for the total D3-D7 system. 

If we concentrate on the £± = 0, £3 = £ case of the P-term model, we know that during 
the inflationary phase (Coulomb phase), supersymmetry should be broken by the FI term 
in the D-term. On the brane side, we know from kappa symmetry that non-selfdual fluxes 
on the D7 gauge field in the K3 directions break supersymmetry with the same mass 
splittings p. Thus one finds that the {T — *J 7 ) 2 term in eq. (|2.14j) is actually the g 2 £ 2 /2 
term of the P-term model. Namely, 

1 f n 2 P 2 

iT 12 _L / T- A *T~ — — . (2.17) 
8g 2 7 Jks 2 v ; 
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The R~ 12 dependence is the same as was found in ||. Note that the last term of ( 2,14| ), 
when expanded, actually contains some of the terms arising from the superpotential in the 
P-term Lagrangian. 

The above analysis more or less summarizes the precise dictionary between the brane 
action and the P-term action of [17|. However, as we discussed earlier, the analysis is 



only for a compact K3 and as such would require a more detailed structure when we try to 
compactify the ]R 2 /Z2 to a T 2 /7<2- In the next section we will carry out the compactification 
of the directions orthogonal to the 7-branes in the full F-theory picture which should also 
capture all of the non-perturbative corrections to the system. Furthermore we will also 
generalize the above configuration to incorporate enhanced global and local symmetries by 
considering stacks of D7 and D3 branes on top of each other. 



3. Generalized D3/D7 model with additional global or local symmetries 

In the previous section, we saw how the Lagrangian for a system of a single D3 and a 
single D7 can be derived using the world volume dynamics of the individual branes. We 
also discussed the additional contributions to the action when hypermultiplets are added 
to the system. The Lagrangian that we get is precisely the expected one, studied earlier in 

All the discussion we had in Section 2, however, was only for non- compact R 2 /Z 2 , 
i.e the orthogonal directions to the D7 (and also the D3) branes were still taken to be 
non-compact (the K3 directions, however, were already assumed compact in Section 2). 
As we know, the non-compact cases are only an approximation to the compact case studied 
earlier in § . Going to the compact case will mean to consider the orthogonal space as a P 1 
or T 2 /l2 (where I2 is the orbifold action that reverses the two direction of the T 2 torus) on 
which both the D7 and the D3 branes are points 7 . Of course, charge conservation requires 
that we should have no global seven brane charges floating around, which, in terms of the 
underlying F-theory construction requires us to view the system as a Weierstrass equation 

@ 

y 2 = x 3 + x F(z) + G(z) (3.1) 

with z as the coordinate of the P . The axion-dilaton of type IIB forms the modular pa- 
rameter r of the torus that is determined by the degree eight and degree twelve polynomials 
F and G respectively using the J function. All these have been defined earlier in Q so we 
will not discuss them further here. The global monodromy (which determines the global 
seven brane charges) is zero and the D3 brane charge is cancelled by switching on NSNS 
and RR three-forms satisfying the equation 

2i / ImrwiAw 2 = 23, (3.2) 
JM 

7 Recall that the compact P 1 appears from the F-theory picture that we discussed in Q. In terms of the 
F-theory four-fold, this is simply a KS x K3 compactification with fluxes. 
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where we have chosen a single probe D3 brane on the manifold Ai = K3 x T 2 jXi- The 
relation between lo\ and L02 to the three forms H^s and Hrr is given in Q. From the 
F-theory point of view, the underlying four-fold has an Euler number of 576. We are also 
assuming that the D3 brane is located at a point z\ on the P , whose coordinate is z. 

In the above relation ( |3.2[ ), if uj\ and 0J2 are arbitrary, then in general supersymmetry 
is broken. This would imply that the D3 brane is moving towards the D7 branes. In Q 
we called this phase the Coulomb phase of the hybrid inflation. In the situation where 
we have isolated one D7 from the bunch of seven-branes, the D3 will eventually fall into 
the D7 brane as an instanton. Because of the presence of FI terms, i.e, because of the 
gauge field J-, this instanton will actually be a non-commutative instanton. The condition 
for this final state to preserve supersymmetry can be described as follows. If one isolates 



from the F-theory curve (3.1) one of the singularities of F and G, at which the torus 
fiber (parameterized by a complex coordinate £) degenerates once, and if one assumes the 
existence of a single normalizable (1,1) form Q in the neighborhood of this singularity, then 
supersymmetry will eventually be restored if the four-form 

F A Q - u x A d( - u 2 A d( (3.3) 



is self-dual. Note that the self-duality of the four-form ( |3.3[) now does not guarantee the 
self-duality of the seven brane gauge fluxes, i.e, we have in general J- 7^ *J~, where the 
Hodge star is with respect to the underlying four-dimensional K3 base 8 . We have also 
identified the one forms 

^ dY) = \-m^r^ m^7) (3 ' 4) 

as determining the T 2 fiber of our F theory construction. For more details, the reader may 



want to consult [43|. In our earlier paper [0], this final stage, when the probe D3 brane 
falls into the D7 brane as a non-commutative instanton, was called the Higgs phase of the 
hybrid inflation. To compare this situation to the case without any incoming D3 brane, 
observe that in that situation supersymmetry would be restored when 0J2 = (the Hodge 
star now being with respect to the six dimensional compact manifold). Of course, for our 
case, since we will always have the probe D3 brane in the picture, this situation will never 



arise. Further details on this has appeared in [44]. 

For the compact case that we discussed earlier in [|j the metric for the system can be 
easily written down taking into account all the effects of the backreactions of the branes 
and fluxes on geometry. Following [44], we write the metric for the system as: 

ds 2 = A _1 (isoi23 + A(dspi + ds 2 K3 ) , (3.5) 

where A is the warp factor. The metric looks similar to the metric for a single D3 brane at 
a point on the compact manifold A4, because the fluxes and branes simulate the situation 



8 The fact that there exists a solution with self-dual four form (3.3) is highly non-generic and has been 
discussed earlier in MM. 
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of having D3 branes at a point on A4. The solution to the warp factor is a little involved, 
because we now require the situation where we have isolated (a) one D7 from the whole 
bunch of seven branes, or (see the next subsection) (b) two D7 branes from the bunch of 
seven branes. The second case (b) is actually what we are mostly interested in this paper, 
because, as discussed in earlier sections, this will give rise to two hypermultiplets on the 
D3 brane and also global symmetries in the system. Before we go into the issue of the 
warp factor, let us therefore see how global symmetries arise in the compact model from 
F-theory. 



3.1 Global symmetry from two D7 branes 

There are two issues that arise when we want to incorporate global symmetries to the 
effective world volume theory on the D3 brane. As we will generate these global symmetries 
from the gauge symmetries on coincident D7 branes, the first issue is the decoupling of D7 
brane gauge dynamics from the D3 brane theory. However, what would suffice for our case 
will be to make the D7 gauge coupling constant sufficiently small compared to the D3 brane 
gauge coupling. This is possible for our case because the effective D7 coupling is suppressed 
by the volume of the underlying compact K3 manifold (see eq.(2.7)). Therefore when we 



determine the covariant derivatives on the D3 brane world volume, we can effectively 
neglect the interactions due to the D7 gauge fields, because the K3 manifold is assumed 
to be fixed at a large volume. On the other hand, the kinetic term of the D7 brane will 



survive, and, bearing in mind the correspondence ( 2.17 ), we can use this to study the FI 
and other relevant terms for our system. 

Secondly, as the existence of the global symmetry is related to the existence of a gauge 
symmetry on the D7 branes, this should appear from the Weierstrass equation that we 
discussed in (jjP|), Let us see how this works in more detail. First, observe that this 
situation is not the constant coupling scenario of 4B, 46]. What we require is that the 
discriminant D have a double zero at a point z\ £ 



si. 



D = 4F 6 + 27G = (z — z\) H{z) = , (3.6) 

where H{z) is a generic polynomial of degree 22, which is regular at z\. This would imply 
that, for a contour C(z±) surrounding the point zi, we will have: 

air = 2. (3.7) 

C(zi) 



The requirement on the discriminant D in (eq. |3.6| ) leaves the freedom to choose the 
functions F and G. Although there are possibilities of various choices that allow a double 
zero, one particular set that will be consistent for our scenario would be: 

F(z) = (z- Zl ) F(z), G(z) = (z-z 1 )G(z), (3.8) 

where F and G are polynomials of degree seven and eleven respectively. This means that 
the F-theory curve near the point z = z\ will be: 

y 2 = x 3 + axz + 0z, (3.9) 
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where z is now z = z — z\ and we have absorbed F and G in the definition of a and (3 



respectively. From the above relations and from the parametrizations used in [47|, we see 
a local A\ singularity. 9 This local singularity will then appear on the D3 brane as a global 
SU(2) symmetry. The fact that there could be an SU(2) symmetry that doesn't lie on the 
constant coupling moduli space of F-theory was also observed in |^(|. In fact, the coupling 
r at the point z = Z\ can be seen to be 

2 

r ~ — ln(z - zi). (3.10) 

ZTTl 

Observe now that as z — ► zi, r — * ioo which corresponds to weak coupling, and therefore 
quantum corrections will not modify this behavior. Thus this configuration will survive 
perturbative as well as non-perturbative corrections. 

There are other interesting issues here. As we observed above, the global symmetry is 
governed by the F-theory singularities. Therefore, if we choose the singularities such that 
the symmetry is now E n then we can have theories that have exceptional global symmetries, 
perhaps giving rise to exceptional semi-local strings. For example, a simple way to get £"8 
as a global symmetry is to allow the following behavior of the discriminant D: 

D = ( Z - Zl ) 10 H(z) , (3.11) 

where H(z) is now a polynomial of degree 14. As discussed in [filjfl , this can be realized at a 
constant coupling moduli space of F-theory by keeping ten seven branes at one point. The 
coupling near the system is strong and therefore not realized as a simple orientifold model. 
We will however not go along those interesting directions in this paper. The possibility 
of the existence of semi-local strings in theories with exceptional global symmetries were 
hinted at earlier in J51] ], although no concrete realization of this has yet appeared in the 
literature. For the D3/D7 system, we see that such a scenario may be possible to realize, 



as most of the global symmetries come from the Weierstrass equation in F-theory [46]. 
More details will appear elsewhere. 

Before moving ahead, let us summarize the content of this section in terms of an 
effective gauge theory: In this section we have described a D-brane picture with a local 
U(l) gauge symmetry and a global SU{2) symmetry. This should now be compared with 
the discussion that we had in section 2 for a U (1) gauge theory and a single hypermultiplet. 
In terms of an effective 4-dimensional gauge model, we will have the following potential: 



V = 2g 2 



2 g 2 
+ — 

2 



|$+| 2 + |<l + | 2 - |$_| 2 - |l>_| 2 -£' 2 



+ 



2 g 2 \s\ 2 i$+r + i$_r + + . (3.12) 

where (<5+, form the second hypermultiplet. We are also making a choice = and 
£3 = £ 0, with regard to eq. ( p. 15 ). This corresponds to a choice of the canonical basis 



9 This curve was also observed in [H in Af — 2 gauge theory with two flavors. This is basically related 
to the Argyres-Douglas point |]49|. 
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Figure 1: The hypermultiplets from two D3 branes and one D7 brane. 



for the 2-form T when only and J-%$ are non- vanishing, so that (.T^n) 2 is given by 
(•^67 ~~ J~s$) 2 . In terms of the brane language used in section 2, the above potential can be 
easily derived from the D7 brane and the D3 brane actions simply by taking traces over 
the adjoint representations of the gauge fields and by neglecting the SU(2) D-terms. 

3.2 Enhancement of local symmetries 

So far, we have been taking one D3 brane probing a system of D7 branes. There is another 
interesting scenario where we put two D3 branes probing the region near a single D7 
brane in the compact model. Of course, since the coupling of the D3 brane cannot be 
ignored compared to the coupling of the D7 brane, we cannot regard this as enhancing 
the global symmetry. The global symmetry remains U(l), but the gauge symmetry will 
now increase. For a single D3 probe, when it is very near the D7 brane and 07 plane, 
we know that the gauge symmetry is SU{2) broken to U(l) at all points on the moduli 
space 10 . This SU(2) symmetry is actually USp(2) and therefore doubling the number of 
D3 branes will make the symmetry USp(4), 11 and, more generally, USp(2k) for k D3 
branes. Thus, classically the symmetry is USp(4) x U(l). The former being a local gauge 
symmetry and the latter global. The dynamics of this system under full perturbative and 



non-perturbative corrections has been discussed in []53 1 . 

To see the hypermultiplets in this setup we can use fundamental strings stretching 
between the branes. In Figure 1, we see that there are two different possibilities. In one 
of the possibilities, depicted as a string stretched directly between two D3 branes and the 
D7 brane one obtains two hypermultiplets as desired. For notational purposes we will call 
this state 1. But there are other possibilities for the compact case where, due to global 
seven brane charge conservation, we require definite numbers of orientifold planes and seven 
branes. Now we can have a string that goes all the way across an orientifold plane and 
comes back to the D7 brane (we call this state 2 in the figure). Under non-perturbative 



10 Classically the gauge symmetry is SU(2) at the point where Higgs vev is zero. Quantum mechanically 
there is never an SU (2) enhancement J5^] . 

11 Again, far away from the orientifold plane and D7 brane the gauge symmetry is SU(2), similar to the 
U(l) case for the single D3 brane. 
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corrections, the 07 planes will in general split into pairs of seven branes which we denote 
as (0,1) and (2,1) seven branes, the numbers being the charges measured with respect to 
axion and dilaton, respectively (in this notation, a D7 brane will be denoted as (1,0)). The 
splitting of the orientifold plane into seven branes create branch cuts on the P . Therefore 
the strings reaching the D7 brane undergo changes in their charge distribution pQ] . If we 

call the strings leaving from D3 as f ^ ^ , then the condition for them to reach the D7 

brane will be 

M { 2,l]-M m .(^ = ^y (3.13) 

where Mtp^-i are the monodromy matrices whose explicit form can be taken from p0| . One 
can check that this condition is met for our case. The masses of the strings stretched 
between the D3 branes and D7 brane can be calculated and will in general be much 
heavier than state 1. Some of these states may even be non-BPS. Since all these states are 
much heavier than the state 1 we can effectively ignore those states for our purposes and 
concentrate only on the state 1. 

Before we end this section, let us mention a few words on the cosmology of this model. 
One first notes that the SU(2) coupling of the two D3 branes (at far distance from D7) 
can no longer be much smaller than the U(l) gauge coupling of the D7 brane unless the 
compactification volume is very small. Hence, the SU(2) gauge symmetry can no longer 
be treated as effectively global as we mentioned earlier. D-term inflation for such a more 
complicated gauge group has not yet been studied. Therefore, in this paper, we will focus 
on a D3/D7 brane model with two D7 branes and two hypermultiplets in the effective 
U(l) gauge theory with global SU{2) symmetry. 



3.3 Analysis of the warp factor 

Now that we have seen how global symmetries and hypermultiplets appear in our setup, it 
is time to go back to the issue of warp factors that we briefly mentioned at the beginning of 
this section. Since we have isolated two D7 branes from a system of branes and orientifold 
planes, we have to carefully consider the equation of the warp factor. In the notations that 
we used above, we keep a probe D3 at a point z = z\ on the P . The two D7 branes are 
located at z = z\. The metric ansatz that we used in fl3.5|) involves two things: (a) the 
metric on the sphere P , and (b) the warp factor A (we will not be too concerned with the 
metric of K3 here). Let us first tackle the metric on the P 1 . In terms of the eta- function 



rj we can express the metric near two D7 branes on top of each other as [54]: 



ds^ = g z2 \dz\ 2 ~ Im t(z) \t] 2 (t) (z - zi)~^dz\ 2 (3.14) 

Putting a probe D3 brane and additional fluxes only introduces a warp factor A to the 
system as discussed above in ( |3.5[ ). The equation for the warp factor can be approximated 
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near the D7 branes as: 

*□ A 2 = H NS A Hrr + (4vr 2 a / ) 2 S 2 (z - z x ) 5 4 K3 + 

+ ^r-[tr (TCAft)-tr (FAf)] 5 2 (z- Zl ) 
64 

where the factor of ^ comes from the gravitational couplings on the two D7 branes. In 
writing the above relation we have ignored the effect of other seven branes that are placed 
far away from the two D7 branes. The position of the D3 brane on the K3 is denoted 
symbolically by the delta function <5j^ 3 . 

The solution to the above equation involve various things: (a) the choice of the three 
forms Hns an d Hrr, (b) the choice of the world volume gauge fluxes J-, and (c) the 
curvature tensor 1Z. The full solution to this equation has not been worked out, and 
therefore we will not try to address it here in details as this would lie beyond the scope 
of this paper. We only note that an approximate solution to the warp factor can be given 
in terms of Green's function a la [p^j . Using this one can show that the warp factor 
is typically of order 1 here as the additional factors are suppressed by the volume of the 
compact manifold. Since we take large internal volume (which earlier helped us to decouple 
the seven brane gauge interactions) we can safely ignore the warping effects here. Therefore 
instead of going to a more detailed analysis of the warp factor, we now turn to the analysis 
of semilocal strings in this model. 



(3.15) 



4. Semilocal BPS-strings in supergravity 

Semilocal strings are a special class of vortex solutions that were first found in the extended 
Abelian Higgs model |3l], [3^ ]. In that model, the Higgs sector of the conventional Abelian 
Higgs model (which gives rise to the standard ANO-strings) is doubled so that the theory 
has a global SU{2) g symmetry, in addition to its local U(l)i gauge symmetry. 12 In the 
last few sections we saw in detail how the D3/D7 system realizes easily the required global 
and local symmetries of this model, which is the main requirement for the presence of the 
semilocal strings. However, our model has much more structure, as it is also embedded 
into N = 1 supergravity. In this section, we make use of this extra structure and show how 
the semilocal strings can arise as 1/2-BPS solitons in the extended Abelian Higgs model 
coupled to four-dimensional, M = 1 supergravity. While it is known that semilocal strings 



can satisfy BPS-type energy bounds in the case of critical coupling |3l|, g^, [34}| , the precise 
connection with local supersymmetry has so far not been established in a clear an precise 
way. 13 Our treatment extends the discussion of ref. [ 23 1 for the locally supersymmetric 



ANO-vortex to the semilocal cosmic string. In the analysis below, we will ignore the 
effectively decoupled D7 brane gauge fields completely, dropping their kinetic term. 



12 Simple generalizations to N Higgs fields with a global SU(N) g symmetry appeared in |33|, whereas 
5l| discussed other symmetry groups that could possibly also give rise to semilocal defects. 
13 For related work in three dimensions, see |3ql. 



4.1 The Lagrangian 



In section 2, we derived the essential features of the low energy effective action of a simple 
D3/D7 brane system in a K3 x T 2 /Z2 compactification of type IIB string theory (see eq. 
Q2.15D for the M = 2 supersymmetric non-gravitational part). As sketched in section 3, 
this system can be generalized to describe one D3 brane near two coincident D7 branes 
with gauge group U(l)i and an effectively global symmetry SU(2) g . This global symmetry 
acts by rotating two (rigid) J\f = 2 hypermultiplets, ($+,$_) and ($+,$_), into each 
other 14 . The two hypermultiplets are accompanied by an J\f = 2 vector multiplet (W^,S) 
corresponding to the gauged U(l)i symmetry (we are suppressing the fermions for the 
moment). It can be shown that the neutral and negatively charged scalar fields have to 



vanish on the string solution [57, 58 1. So we will from now on set in this section 

S = = l>_ = (4.1) 



in equation (3.12). This collapses the superpotential, which we will henceforth set to zero. 
For the two remaining, positively charged, scalars and we employ the somewhat 
more convenient notation 

ipi = $+, w = $+ . (4.2) 



Coupling this model to gravity reduces supersymmetry to M = 1, and the total field 
content including the fermions now derives from the M = 1 supergravity multiplet (e^jip^), 
two chiral multiplets {xii^i) (i = 1>2) and a vector multiplet (W^,X). Clearly, the two 
chiral multiplets carry an electric charge +1 under the local U(l)i symmetry gauged by 
Wfj, and transform as a doublet under the global SU(2) g symmetry. In the following, we 
employ the supersymmetry conventions of ref. |l6|, 56 1. In particular, we use (p l to denote 
the complex conjugate of ipi and indicate the chirality of their superpartners by Xi (left 
handed) and x % (right handed). 

Taking the Kahler potential and the gauge kinetic function to be minimal, the bosonic 
part of the Lagrangian then reads 



e-^bos = -\M%R - D^i £>V - IF^F^ - V D , (4.3) 
where the D-term potential is defined by 

V D = \d 2 D =g(t-<p i <p i ) (4.4) 

and 

= d^W v - d v W» , D^i = (0„ - igW^tpi (4.5) 

14 SU(2) g should not be confused with the R-symmetry group SU(2)r, which rotates <J?+ into $_ and 
3>+ into 
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denote the usual Abelian field strength and the U(l)i covariant derivative. The transfor- 
mation rules of the fermions in a purely bosonic background are 

S^L = (fy + |^ a6 (e) 7o6 + |L4j) e L , 
5xi = z(0-ig W)W£R, 

SX = \^ u F^ v e + ii 75J De . (4.6) 
Here the (composite) gravitino U{\) connection, Aj^, is given by 

A * = 2M 2 } ~ + W W " D 

4.2 Semilocal vortices with unbroken supersymmetry 

Our goal here is to show how the above action can give rise to string-like defect solution 
that preserve one half of the original J\f = 1 supersymmetry. For simplicity, we restrict 
ourselves to non-singular configurations that are static and cylindrically symmetric, that 
is, we assume (1 + l)-dimensional Poincare symmetry on the 'world sheet' of the string as 
well as a rotational symmetry with orbits orthogonal to that world sheet. Using (t,z) as 
the world sheet coordinates, we align the rotation axis of the string with the z-direction. 
The orthogonal distance from the rotation axis is denoted by r and rotations about the 
axis are parameterized by an angular coordinate 9. The space-time metric can be brought 
to the form 



ds 2 = -dt 2 + dz 2 + dr 2 + C 2 (r)d8 2 . (4. 
We use the vierbein with e 1 = dr and e 2 = C(r)dd, which leads to the spin connection 



0, u e 12 = -C'(r). (4.9) 



In order for the string to have finite energy per unit length, the potential and kinetic 
energy densities have to vanish sufficiently fast at large distances r from the string axis. As 
for the potential energy, this means that the scalar fields (fi have to approach the vacuum 
manifold S 3 = {(fi : \<fi\ 2 + |^2| 2 = £} of the potential V D as r —> oo. A cosmic string in 
this model thus defines a mapping of the spatial circle at infinite r into the vacuum 
manifold S* 3 . 

If this was the only condition, one would conclude that there could be no topologically 
non-trivial vortices, as 7Ti(5 3 )={l}, i.e., any non-trivial mapping of could be deformed 
into the trivial map (i.e., the vacuum) at no cost of potential energy. 



As was pointed out in [31, 32], however, the requirement of finite kinetic energy per 
unit length is of equal importance, as it implies that D^ifi also vanish sufficiently fast at 
large r. This further restricts (<pi,<f2) to approach a U(l)i gauge orbit (i.e., a great circle 
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Sf) on S 3 as r goes to infinity. The Higgs field (<pi,<p2) at r — ► oo thus really defines a 
map of 5^ into a gauge orbit S 1 of U(l)i on the vacuum manifold S 3 . Any such map is 
characterized by a winding number n G TTi(Sf) = Z and an associated gauge connection 
that compensates the gradient energy from the angular dependence of the Higgs field. 
Stokes's theorem and the requirement that Dnfi — > as r — > oo then imply the presence 
of a magnetic flux through the (r, 0) plane: 

W d d6 = 2irn. (4.10) 

As there are no gauge orbits perpendicular to the U(l)i gauge orbits (because SU(2) g 
is not gauged), any collective motion off such a U(l)i orbit (i.e., any unwinding of the 
string at infinity) would cost an infinite amount of gradient energy 15 which could not be 
compensated for by an appropriate gauge connection. Thus, it is really the fundamental 
group of a gauge orbit Sf on the vacuum manifold that matters for the stability against 
unwinding and not the naive fundamental group of the vacuum manifold S 3 itself 32]. 



It is important to note, however, that the stability against unwinding (i.e., the con- 
servation of the total magnetic flux as a topological invariant) does not necessarily mean 
that such defects are stable against other types of deformations. In particular, the mag- 
netic flux need not be confined or localized in a small or thin region of spacetime, but may 
spread out to fill larger regions (while still being globally conserved). As was pointed out in 
|32] , |33| , |35[| , topological defects may be unstable against such 'spreading' if the full vacuum 
manifold is a non-trivial fiber bundle over a base manifold with the gauge orbit as the 
fiber such that the relevant homotopy group of the total bundle becomes trivial. To trans- 
late this to our situation, one notes that S 3 is a nontrivial bundle over CP 1 = S 2 , and 
tti(S 3 ) = {1}, whereas Tri(Sf) = Z. 16 As we will see shortly, in our model this phenomenon 
manifests itself in the form of zero modes that describe a possible fattening and warping 
of the magnetic flux tube. This ability to spread out the flux is one of the properties that 
make these semilocal strings more appealing in a cosmological context, see the discussion 
in section 5. 

Let us now describe the semilocal string solutions of our model in more detail. The 
above-mentioned asymptotic behavior for r — > oo and the cylindrical symmetry imply that 
$(r, 6) = ((pi(r,6),p2(r,0)) approaches a field configuration of the form 

$(r,0) — ><S> e in6 , r — > oo (4.11) 

where = £• The global SU(2) g symmetry allows us to align $o with the if \ direction, 

and the most general Higgs field configuration consistent with cylindrical symmetry is 
then 17 

In a finite space, this gradient energy barrier would be finite. 

16 In |Q such a nontrivial fibering with a collapse of the relevant homotopy group over the total bundle 
has been suggested as the defining property of a semilocal defect. In the defining property is chosen 
to be the ability to 'spread' the magnetic flux. 

17 Linear combinations such as ~^2 m hS m ' (r)e lm6 would break the assumed rotational symmetry, as quan- 
tities such as the potential V D would in general no longer be ^-invariant. 
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(4.12) 



/ Mr,0)\ = [ f(r)e ine \ 
ytp 2 (r,e)J \h{r)e imd )' 

whereas the gauge potential takes the form 

gW u dx^ = na(r) d9 -» F = \F UV dx» dx u = ^^dr d9 = ^^rf^eV . (4.13) 

9 gC{r) 

Here, f(r), h(r) and a(r) are the profile functions of the Higgs fields and the vector 
field, whereas m is an a priori arbitrary integer. 

The desired asymptotic behavior at r — > oo and the assumed non-singular behavior at 
the origin impose the following boundary conditions on the functions f(r), h(r), a(r) and 
the metric function C(r): 

C(0) = 0, C"(0) = 1 
/(0) = 0, /(r) — sfi as r —> oo 

h(0) = ho5 m o, h(r) — > as r — > cxd 

a(0) = 0, a(r) — > 1 as r — > oo (4.14) 

where we have indicated that h(0) might be non-zero for m = 0, with ho being some 
constant. 

Our goal now is to find profile functions f(r), h(r), a(r), C(r), which describe a 
solution of the Lagrangian that preserves 1 /2 of the original supersymmetry. We thus have 
to find an appropriate Killing spinor e(r, 9) such that Sip^ = 5\i = SX = can be satisfied 
for a certain form of the profile functions. For simplicity, we will from now on assume 

n>0. (4.15) 

In order to cancel the magnetic flux contribution to 5X, we choose the following projector 
for the Killing spinor 

7 12 e = -i 75 e • (4.16) 



The BPS equations that follow from the vanishing of the transformation rules ( |4.6| ) for 
the chiral fields \% an d for the gaugino A are then 

(C(r)d r + i(d e -igW e j)<pi = 0, 

F 12 -D = 0, (4.17) 



which leads to 



C(r)f'(r)=n[l-a(r)]f(r), (4.18) 
C(r)h'(r) = [m-na(r)]h(r), (4.19) 

a'(r) = °-^l[H-f\r)-h 2 (r)]. (4.20) 
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These three equations are not completely independent, as solving eqs. ( 4.1§| ) and ( 4.1E| ) 
for a(r) and equating both expressions yields the relation 18 



d r In 



h(r)\ m — n 



C(r) ■ 



(4.21) 



We make further use of this relation in Section 2.3 to derive an upper bound on the possible 
integers m. 

It remains to determine the gravitino BPS equation, for which one needs the composite 
connection : 



A? = 0, Af 



M 2 



[nf 2 (r) + mh 2 {r) + na{r) (f - f(r) - h 2 (r))] . (4.22) 



The gravitino BPS equation is now 

d r e = , [d e + \iC\r) + ±L4|] e L {6) = . 



(4.23) 



The 9 dependence of e follows from the limit r — > 0, where Aq goes to zero and C'{r) 
approaches 1 so that ( |4.23| ) implies 



e L (9) = e-2 ie e 0L 



(4.24) 



where eoL is the left handed part of a constant spinor satisfying the projection equa- 
tion ( f4,16| ). Thus, the gravitino equation dip^ = boils down to 



1 _ C'{r) = A 



B 



(4.25) 



with Af as in ( ^22| ). 



4.3 Limiting cases 



While no analytic solution to the BPS equations ( 4.18| ) - ( 4.20| ), (4.25) is known, a great 
deal of physical information can already be obtained by studying the limiting behavior at 
r = and r — > oo. 

We will first show how the limiting behavior constrains the possible integers m in 
(p2(r,0). In the limit where r approaches zero, C(r) = r + 0(r 2 ) (see eqs. ( 4.14|) ). As a(r) 



must be at least of linear order in r, the a-term in the BPS equations (|4.1g| ) and ( |4.1£| ) 
can be neglected to lowest order in r. To lowest non-trivial order, these two equations thus 
read 



rf'(r) = nf(r) 
rh (r) = mh(r) 



(4.26) 
(4.27) 



18 In the special case n = 1 and m = 0, the flat space limit of (4.21) reduces to the expression h(r) = 
(w/r)f(r) used in [Q. Analogues in conformally flat worldsheet coordinates also appear in J34|. 
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whence 

f(r) = f n r n + 0(r n+1 ), h(r) = h m r m + 0(r m+1 ) (4.28) 

with some non- vanishing coefficients /„, h m (h m might be zero for m = 0). Finiteness of 
h(r) at r = then implies m > 0. 

To derive also an upper limit on m, consider eq. ( [4.21 ). As h(r)/f(r) — > for r — > oo, 



the left hand side of ( [4.21 ) has to be negative, which, because of the positivity of C( 



implies that m < n. We have thus obtained the bounds 

< m < n. (4.29) 

Let us now take a closer look at the limiting behavior at large r. As r approaches 
infinity, (/(r), h(r),a(r)) approach their asymptotic values 0, 1). It is easy to see, that 
the four BPS-equations are consistent with this behavior provided that asymptotically 

C(r) ~ r ( 1_ i|)- (430) 

This shows that asymptotically the string creates a locally-flat conical metric with an 
angular deficit angle which is due to the constant FI term £: 



ds z = -dt 2 + dz z + dr 2 + r 2 ( 1 - d6 2 . (4.31) 



Notice also that in the limit r — ► oo the full supersymmetry is restored because = 0, 
D = 0, d r ifi = Dgtpi = and Rf_ iu ab = which corresponds to the enhancement of 
supersymmetry away from the core of the string. This is the same result as for the standard 
ANO vortex |2j|, however, it should be pointed out that the rate at which the asymptotic 
values of the profile functions f(r), h(r) and a(r) are obtained is in general no longer 
described by an exponential fall-off, but can change into a power law behavior depending 



on the values of certain moduli parameters of the semilocal string [32| . To understand 
this better, let us now turn to the behavior at small r. 

Expanding the profile functions as power series of the form f(r) = /o + fir + f2f 2 + . . . 
and similarly for h(r), a(r) and C(r), we know already that 

f(r) = f n r n + 0(r n+1 ) (4.32) 
h(r) = h m r m + 0(r m+1 ) (4.33) 
C(r) = r + 0(r 2 ) (4.34) 

In order to determine the lowest non-trivial power of a(r), we expand (|4.20|) to lowest order 
and obtain 

a 2 

a{r) =a 2 r 2 + 0{r 3 ), with a 2 = (4.35) 

In 

In a similar fashion one can now determine the higher coefficients order by order and 
arrives at recursion relations that determine all coefficients in terms of £ and the lowest 
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non-trivial coefficients f n and h m of the Higgs fields. It is easy to see that the recursion 
relations are such that f(r) contains only odd powers when n is odd and only even powers 
when n is even. Similarly, the function h{r) is odd/even if m is odd/even. a(r), by contrast 
is always even, whereas C(r) only contains odd powers. 

This power series expansion leaves the lowest Higgs coefficients f n and h m undeter- 
mined. In fact, the coefficient h m turns out to be a undetermined parameter of the solution, 
whereas f n can be fixed numerically by matching to the required asymptotic behavior at 
infinity. The numerical value of /„ depends upon the parameter h m . 

The parameters h m can actually be taken to be complex, as one has the freedom 
to shift the phase of ip2 = h(r)e imd by a constant shift and absorb that phase into h m . 
If one breaks the rotational invariance by either allowing linear combinations (f2(r,8) = 
Xym=o h( m ~) (r)e tm ® and/or by splitting an n- vortex into (up to n) vortices with smaller 
winding numbers and different positions in the (r, 6) plane, one expects supersymmetry 
to be preserved and finally arrives at a 2n-complex-dimensional moduli space of semilocal 
vortices [34]. 



The physical interpretation of the moduli h m can most easily be illustrated for the 
parameter ho in the case m = 0. In that case, eqs. ( 4.35| ) and ( 4.1 3| ) yield the magnetic 
field 

F 12 (0) = <?(£- (4.36) 
at the origin. In the limit ho — > this approaches the magnetic field in the core of the 



ANO-string [23]. This is not surprising, as ho = implies h(r) = in the case of m = 
via the recursion relations that follow from the BPS equations. In that case, the BPS 
equations reduce to those of the ANO-string. 

On the other hand, increasing /i from to decreases the magnetic field in the 
core. As the total flux only depends on the winding number n, and thus is unaffected by 
changes in ho, this means that the flux tube must become flatter and wider. Thus, ho is a 
measure for the width of the magnetic flux tube and encodes the possibility for semilocal 
strings to spread out their flux into larger space-time regions |p2| . 

Another interesting quantity is the spacetime curvature R = 2C" /C in the center of 
the string. At r = 0, this quantity reduces to R(0) = I2C3. C3 can be easily obtained from 
the first non-trivial power of eq. ( 4.25| ), which leads to 



m = _ i iA\ _ 2 (9iL_3l\ (437) 



M P J \ M P 

We finally note that for high n and m, the functions f(r) and h{r) are zero at least up 
to order r m (remembering m < n). Up to terms of order r 2m , the BPS equations ( 4.2C] ) 
and ( f4.25| ) then simplify to 

a'(r) = —C(r) , 1 - C'(r) - nMp 2 £a(r) = , (4.38) 
n 
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which, in view of the boundary conditions at r = 0, can be solved as 

a(r) 
C(r) 



M 2 P ( gt . 



M P g£ 
■ sm( r) 



Up to terms of order r 2m+1 , the metric at small r is then: 



ds 2 



-dt 2 + dz 2 + dr 2 + 



whereas, up to terms of order r m for the Higgs and r + for the gauge fields, one obtains 



M P gi 
— — sin — — r) 

Zq y M P 1 

2m+2 



de 2 



(4.39) 



(4.40) 



F re = M P sm(^-r). (4.41) 



5. Cosmology and semilocal strings 



5.1 Inflation 

The semilocal strings embedded into supergravity were shown in the previous section to 
have an unbroken supersymmetry for a continuous range of the parameter ho, which ef- 
fectively measures the width of the flux tube and is not defined by the parameters in the 
action. Now we would like to study the properties of the semilocal strings in applications 
to cosmology. 

The generalized D-term model after doubling the number of hypermultiplets corre- 
sponding to doubling the number of D7 branes in the D3/D7 model has the following 
potential: 



V 



g_ 

2 

2 1 1 



(f>+4>- - 



+ 



9 



M 2 + 



+ 



+ i0+r -10- 1 - 10- 



+ 



2^ 



+ 



(5.1) 



Here all scalar fields have canonical kinetic terms, such as ^(dcj)) 2 . They are related to 
the chiral fields in the action (5.1) as follows: 0+ = \f2<& + ,s = y/2S etc. The tree-level 
potential ( |5.1| ) has a supersymmetric minimum (a valley with V = 0): 



s = 

I i |2 







+ 



2£ 



(5.2) 



The U{\) gauge symmetry is broken in the minimum. In the absence of strings, the fields 
can take any values along this valley. When strings are present, the fields tend to drift 
towards the state with 0_ = 0_ = 0, because <^>_ = 0_ = for the BPS (lowest energy 
state) solutions [57, 55]. 
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The potential has a second flat direction, 



0. Inflation occurs when the field 



s rolls along this flat direction. At the classical level, the scalar potential (|Tl]) equals 4>-£ 2 
for all values of s at <f>± = <f)± = 0; in order to study inflationary dynamics one should take 
into account one-loop corrections to the effective potential. 

We will consider the models with g 2 ^ 10 -4 , because it is more difficult to motivate 
much smaller coupling constant g 2 in string theory. As we will see, inflation in the model 
with g 2 p> 10~ 4 ends at s > s c = ^/t;. In this regime, the one-loop corrected potential of 
the field s in the model with C hypermultiplets of equal g 2 is given by [0, 30] 



a 2 



1 + 



8vr 2 



In % 

si 



To study inflation, one should use the Friedmann equation 



H 



a ^6 



where a(t) is a scale factor of the universe. This leads to inflation 



a(t) = a(0) exp 



During the slow-roll regime, the field s obeys equation 3Hs 

CgHt 



-V'{s) 



s\t) = s 2 (0) 



2^/6vr 2 



(5.3) 
(5.4) 

(5.5) 
which gives 

(5.6) 



Using equations (|5.5[), (|5.6|) one can find the value of the field s N such that the universe 



inflates e N times when the field rolls from s N to the bifurcation point s = s c 



2vr 2 



Cg 2 N 
~2^~ 



(5.7) 



Co 2 

In this paper we will concentrate on the case > £. Density perturbations on the scale 



of the present cosmological horizon have been produced at s ~ sn s c with iV ~ 60, and 
their amplitude is proportional to ^ttt- at that time ISO]. One can find £ using the COBE 



normalization for inflationary perturbations of the metric on the horizon scale [27]: 

1 ^3/2 



This yields 



y3/2 



5V3vr V 
2y/2ir 2 t 



1.9 x 10" 



2tt£i 



5.2 x 10~ 4 . 



V Cg 

In what follows, we will concentrate on the case with two hypermultiplets and take ./V 
In this case 



5.2 x 10" 4 
2^ 



y/UjN ~ 1.5 x 10" 



(5.8) 

(5.9) 
= 60. 

(5.10) 
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This yields the scale of spontaneous symmetry breaking 

y/l ~ 9.3 x 10 15 GeV. (5.11) 

This model has a nearly flat spectrum of density perturbations, n ~ 0.98, and a very 
small amplitude of gravitational waves. Both of these features are in a good agreement 
with the observational data. 



5.2 Semilocal strings in cosmology 

After the field s passes the bifurcation point s c = the fields <f>±,<f>± =0 start rolling to 
the minimum of the effective potential at \4>+\ 2 + |</>+| 2 — |0-| 2 — \<ft-\ 2 = 2£. This leads 
to production of cosmic strings. One should check whether the semilocal cosmic strings 
produced during this process can affect the amplitude of large-scale density perturbations. 

Cosmic strings lead to perturbations of the metric with flat spectrum if in the course of 
the evolution of the universe there always remains at least one string of length comparable 



to the size of the particle horizon ~ ct [61]. We will call such strings infinite. Thus one 
needs to check whether infinite semilocal strings can form and survive long after the end 
of inflation. 

This issue is rather nontrivial. Even though stability of semilocal strings is not pro- 
tected by topological considerations, one could argue that infinitely long straight strings 
are protected because they are BPS solutions, so their energy cannot get any lower. On the 
other hand, the fact that the thin and thick BPS strings have the same tension indicates 
that strings can be unstable if it is possible to excite a zero mode which would lead to the 
growth of their thickness; see e.g. p2| , |33| , |34| , However, it would require infinite action to 
excite such a zero mode in an infinitely large universe [Q|. Also, the growth of thickness of 
a single infinite string occurs with the speed smaller than the speed of light, and the string 
tension does not change in this process. Meanwhile (however strange this might seem), the 
diameter of the observable part of the universe grows with the speed 4 times greater than 
the speed of light in a universe dominated by ultrarelativistic particles, and with the speed 
6 times greater than the speed of light in the cold dark matter dominated universe |60|[ . 
Thus, on the horizon scale, each separate infinite semilocal string will always look thin, and 
its expansion will not lead to the change of its tension. This means that the spreading of 
each individual BPS string would not change the amplitude of density perturbations which 
it produces. 

However, the strings created at the end of inflation initially are not straight, infinite 
BPS strings. In contrast to the case of the standard Nielsen-Olesen cosmic strings, infinite 
semilocal strings do not appear at the moment of the phase transition. Typically, soon 
after the phase transition, strings look like a cloud of string segments of a size comparable 
to the correlation length at the moment of the cosmological phase transition. These string 
segments have global monopoles at their ends. Each global monopole attracts a nearby 
global anti-monopole with a force that does not depend on the distance between them. 
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If the monopole and the antimonopole closest to it belong to the ends of two different 
strings, they attract each other, and the two strings tend to merge, forming one longer 
string. However, if a string is shorter than the distance between it and another string, the 
two monopoles on its ends are attracted to each other stronger than to other monopoles. 
As a result, either such short string segment shrinks and disappears due to its tension and 
monopole interactions, or in the beginning it forms a small loop, and then the loop shrinks 
and disappears. 

Infinite semilocal strings responsible for the large scale density perturbations can ap- 
pear only because of merging of short string segments. The probability of this process is 
expected to be small; it is not reinforced by topological considerations. In addition, the 
possibility of growth of thickness of the near-BPS string segments eventually may lead to 
their overlapping with other string segments, which may result in their de-localization and 
effective disappearance. (Note that it does not cost infinite energy or infinite action to 
excite the zero mode for each segment of a semilocal near-BPS string of a finite size.) A 
numerical investigation of this process was performed in |37|| , with the conclusion that in- 
finite semilocal BPS strings are not formed, and therefore they do not lead to the problem 
of excessively large string-related metric perturbations. 

One should exercise some caution before applying the results of Ref. [ml] to the model 



( |5.1| ). Indeed, numerical simulations performed in Ref. [37] did not fully take into account 
the quantum dynamics of spontaneous symmetry breaking in hybrid inflation. Instead of 
that, the authors simply made some plausible assumptions about the random distribution 
of the scalar field after the phase transition. Then they verified that their results are not 
very sensitive to these assumptions, but in their analysis they did not take into account the 
expansion of the universe. Also, the numerical analysis was performed in |37| for a simpler 
model, which did not contain the fields 4>-,<p-. One could argue |30| that these fields are 
not important because eventually the fields drift towards the state with = = since 
= 0_ = for the BPS solutions 57, 58 1. The short string segments produced after 
inflation initially are not in the BPS state, and the process of disappearance of the fields 
<fi- = (f>~ = could take a very long time but the numerical simulations performed in 



[59] suggest that this time can be rather short. 



We do believe that the results obtained in [37, 59 1 provide a strong evidence that the 
semilocal strings formed in the model fl5,l|) do not pose any cosmological problems, but the 
importance of this subject warrants its full quantum mechanical analysis. This can be done 
using the methods of Ref. |0S| 1, where the process of spontaneous symmetry breaking in 
the simplest versions of hybrid inflation was investigated. These methods allow to perform 
a complete quantum field theory investigation of string formation in the early universe 
]62| ]. Leaving this full analysis for a future investigation, we will discuss some important 
features of the processes which occur after inflation in the model ( |5.1| ) . We will argue that 
the expansion of the universe does not lead to qualitative changes in the processes studied 
in p7f , and that in the first approximation one can indeed ignore the effects related to the 
fields </>_,(/>_. Our analysis will provide a framework for the future numerical investigation 
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of this issue. 



5.3 Symmetry breaking and string formation at the end of inflation 



According to [17, pQ| , the field s at the last stages of inflation rolls down along the valley 
with 4>±,4>± = until it reaches the bifurcation point s c = Then it continues rolling 
along the ridge with (f)±,4>± = until the fields fall down from the ridge and sponta- 
neous symmetry breaking occurs. When the fields 4>±,<j)± acquire nonvanishing (complex) 
expectation values, their initial distribution looks very chaotic, which can lead to string 
formation via the Kibble mechanism. 

Even though the discussion of spontaneous symmetry breaking can be found in ev- 
ery book of quantum field theory, the dynamics of this process was understood only very 
recently, and the results were rather surprising Q. For example, one could expect that 
when the fields fall down from the top of the effective potential to its minimum at <p ~ \/2£, 
they should enter a long stage of oscillations with the amplitude ~ \/2£- However, a de- 
tailed investigation of this issue shows that because of the exponential growth of quantum 
fluctuations in the process of spontaneous symmetry breaking, the main part of the kinetic 
energy of the rolling fields becomes converted into the gradient energy of the inhomoge- 
neous distribution of these fields within a single oscillation. The resulting distribution 
looks like a collection of classical waves of scalar fields moving near the minimum of the 
effective potential, the amplitude of these waves being much smaller than the amplitude 
of spontaneous symmetry breaking y/2£. In other words, the main part of the process of 
symmetry breaking completes within a single oscillation of the field distribution p6|]. 

In order to study this process in the context of the model ( |5.l|) , we need to know the 
speed of rolling of the field s as it passes the bifurcation point s c . This will allow us to find 
how far along the ridge the field s can go until spontaneous symmetry breaking occurs. 
The estimates we are going to make will be very rough, but sufficient for our purposes. 

According to (5.7) and (5.10), for C = 2 and g 1 10~ 4 the slow-roll regime ends 



not at the bifurcation point, but earlier, at s\ ~ - 3> s c . At s < s± the field rolls down 
very fast, reaching the bifurcation point s c within the time smaller than the Hubble time 
H . Therefore during this time the total energy density of the field s is (approximately) 



g A e ln s 2 

bifurcation point. 



- i2 ^-^i In \t for the velocity of the field s at the 



Symmetry breaking happens at s < s c because the fields </>+, 4> + acquire the tachyonic 
mass at s < s c , 



m% = g z (s z - si) < . (5.12) 

This means that, after passing the bifurcation point, the tachyonic mass of the fields 4> + ,(f) + 
grows as follows: 



m\ ss -2 5 2 v^As = g 2 v ^sAt , (5.13) 



where As = s r — s. 
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The process of spontaneous symmetry breaking occurs because of the exponential 
growth of quantum fluctuations of the fields with a tachyonic mass. Ignoring, in the first 
approximation, the time dependence of the tachyonic mass, one can find that quantum 
fluctuations of the fields with m\ < and momentum k < \m+\ grow exponentially: 



~ e 



Occupation numbers of quantum fluctuations grow as e 2 Vl m +l _fc2 *. When the occupation 
numbers become much greater than 1, growing quantum fluctuations can be interpreted as 
a classical scalar field. This means that the growing components of the scalar fields cp + , 4> + 
will look inhomogeneous on the scale Im+j -1 . The average initial amplitude of quantum 
fluctuations with k < \m+\ is O(10 _1 |m + |) |3{J. Because of the exponential character of 
the growth, it takes the field distribution the time ~ |m + | _1 ln before it reaches the 
minimum of the effective potential. In fact, the total duration will be somewhat smaller 
because the curvature of the effective potential grows in the process. 

Note that in our case \m+\ is vanishingly small at the bifurcation point; it becomes 
large only later. Therefore spontaneous symmetry breaking occurs only when the time 
that it takes for the rolling field s to reach some point s < s c becomes comparable to 
|m + (s)| _1 In j^q^fyy- This condition reads 



„2 



, 2/ 7 A] - 2 wvz * 2 9 A e^^ 

To find an approximate solution of this equation, we will take |m+(s)| = |m+(0)| = gy/£ 
under the logarithm. This yields 

As~F(g) g 2 / 3 ^, (5.16) 

where 

Numerical evaluation of this expression shows, for example, that F ~ 0.8, and As ~ 0.2 s c 
for g 2 ~ 10~ 2 . Finally, we find the absolute value of the tachyonic mass of the fields </>+,</>+ 
at the moment when the symmetry breaking occurs: 



2 1 l A F^ 2 (g)g^^l = 0(// 3 ^) • (5-18) 



Note that at this stage the only fluctuations experiencing the tachyonic growth are 
those of the fields (f>+,(f>+; the fields (j)-,(fi- are not generated. On average, the amplitude 
of the complex fields 0+ and should be the same, but in different parts of the universe 
separated by a distance more than m^ 1 the values of these fields as well as their phases can 
be different. Therefore the correlation length which determines the typical initial length of 
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the string segments and the typical initial distance between the cosmic strings is expected 
to be 



Ax = O(3 4/3 70 . (5.19) 

For the range of parameters we are interested in, this correlation length can be two orders 
of magnitude smaller than the horizon size at that time, H ~ g^ 1 £ , which means that 
the corresponding momenta are two orders of magnitude greater than H. 

This implies, in particular, that one can neglect the expansion of the universe at the 
first stage of string formation. This removes one of the concerns related to the applicability 
of the results of Ref. |37| to the investigation of string formation after hybrid inflation. 

The simple tachyonic description of the process of spontaneous symmetry breaking is 
valid only at the first stages of the process. The effective potential of the fields s and </>+,</>+ 
far away from the ridge 4> + ,<p + = is very complicated, so one could expect that the fields 
will move down along a very complicated trajectory. Rather surprisingly, the field s and 
the (quasi)homogeneous component of the fields <^+,(/>+ follow a linear trajectory. Once 
the amplitude of the fields ^>+,<^+ grows up to O(As), they start falling down along the 
straight line 

4>+ = f4>+ , 



1 + f 2 r- 

—±- <f> + + s = , (5.20) 

where / is an arbitrary constant depending on the local value of the ratio between <f> + and 
(f) + when they approach the linear trajectory. A similar regime for the nearly homogeneous 
component of the scalar fields was found in |p3| in a model containing only two time- 
dependent fields, cj) + and s. 

Eventually the fields s, 0+, 0+ roll to the minimum of the potential at s = 0, 4>\+4>\ = 
2£ and start oscillating there. When these fields reach the minimum and continue moving 
forward by inertia, a certain combination of the fields 0_ , 0_ acquires a tachyonic mass, 
and the corresponding fluctuations begin growing exponentially. However, this happens 
only during a very short period of time, of the same order as the inverse tachyonic mass 
of the fields (without the important logarithmic factor In t^t ? which appears in 

the description of the growth of the fluctuations of the fields </>_,0_). As a result, the 
amplitude of the quantum fluctuations of the fields 4>- , 0- does not grow much and remain 
much smaller than the amplitude of the fields 0+ , <^>+ ~ -y/f . 

One should note that this description is a bit oversimplified. Strictly speaking, the 
solution ( |5,20 ), is valid only for homogeneous fields. In our case, the homogeneity is 



broken by tachyonic fluctuations on the scale (<? 4//3 \/f) _1 > but one can expect that, in 
the first approximation, Eq. ( |5.20| ) should still hold in each domain of size (<? 4 / 3 
because the curvature of the effective potential along the main part of this trajectory is 
higher than (ff 4 ^ 3 \/f) 2 - O n a larger scale, the distribution of the fields 0+,^>+ is very 
inhomogeneous. At the time when the fields reach the minimum of the effective potential, 
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the field distribution looks like a collection of classical waves of all scalar fields. These 
waves may move, collide, and produce waves of all other fields, including the fields 0_ , 4>- . 
However, as we already mentioned, a numerical investigation of similar models in [^] shows 
that the amplitude of the waves generated during this process is significantly smaller than 
the amplitude of the fields </>+, 4> + = 0{-\f£). 

This means, in particular, that in the first approximation one can neglect the role of 
the fields (/>_ , 0_ in the process of symmetry breaking and string formation. This is an 
additional argument suggesting that the results of the investigation of string formation in 



Ref. [37] in the model not involving the fields <j>—, (j)- should remain valid for our model as 
well. 

One should emphasize that the process described above is extremely complicated. 
Therefore it would be important to verify our analysis by performing a full numerical 
investigation of symmetry breaking in this model by the methods developed in [39]. We 



hope to return to this issue in a subsequent publication. However, on the basis of our 
analysis we do expect that the main argument contained in [3C] and in our work should 



be valid: By adding a new hypermultipet to the inflationary model discussed above, one 
can avoid the excessive production of metric perturbations related to cosmic strings, while 
preserving the usual inflationary perturbations. 
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